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ABSTRACT 

The age of the universe in the Big Bang model can be calculated from 
three parameters: Bubble's constant, h; the mass density of the universe, 0.^', 
and the cosmological constant, Q./\. Recent observations of the cosmic 
microwave background and six other cosmological measurements reduce the 
uncertainty in these three parameters, yielding an age for the universe of 
13.4 lb 1.6 billion years, which is a billion years younger than other recent age 
estimates. A different standard Big Bang model, which includes cold dark 
matter with a cosmological constant, provides a consistent and absolutely 
time-calibrated evolutionary sequence for the universe. 

In the Big Bang model, the age of the universe, to, is a function of three parameters: h, 
0,m and Q\ (1). The dimensionless Hubble constant, h, tells us how fast the universe is 
expanding. The density of matter in the universe, ^Im, slows the expansion, and the 
cosmological constant, CIa, speeds up the expansion (Fig. 1). 

Until recently, large uncertainties in the measurements of h, and made efforts to 
determine to(^)^m;^A) unreliable. Theoretical preferences were, and still are, often used 
to remedy these observational uncertainties. One assumed the standard model (f^m = !> 
f^A = 0), dating the age of the universe to to = 6.52/ h billion years old (Ga). However, for 
large or even moderate h estimates 0.65), these simplifying assumptions resulted in an 
age crisis in which the universe was younger than our Galaxy (to ~ 10 Ga < tcai ~ 12 
Ga) . These assumptions also resulted in a baryon crisis in which estimates of the amount 
of normal (baryonic) matter in the universe were in conflict (2, 3). 

Evidence in favor of < 1 has become more compelling (4-8), but is still often 
assumed to be zero, not because it is measured to be so, but because models are simpler 
without it. Recent evidence from supernovae (SNe) {4> 5) indicates that JIa > 0. These 
SNe data and other data exclude the standard Einstein-deSitter model {^m = Ij ^^A = 0). 
The cosmic microwave background (CMB), on the other hand, excludes models with low 
Qrri and r^A = (5). With both high and low 0,rn excluded, Q\ cannot be zero. 
Combining CMB measurements with SNe and other data, I (5) have reported 
r^A = 0.62 ± 0.16. [see {10-12) for similar results]. If Oa 0, then estimates of the age of 
the universe in Big Bang models must include Oa- Thus one must use the most general 
form: to = f{^^,^K)/h {13). 



2 



2.0 



K 1.5 



1.0 



0) 

N 0.5 



0.0 



(1.0, 0.0 

(0.3, 0.0 

;o.o, 0.0 

(0.3, 0.7 

'0.1, 0.9 




-20 -15 -10 -5 NOW 
Time [Ga](o,68/h) 



5 



10 



Fig. 1. — The size of the universe, in units of its current size, as a function of time. The 
age of the five models can be read from the x axis as the time between NOW and the 
intersection of the model with the x axis. The main result of this paper, to = 13.4 it 1.6 
Ga, is labeled "to" and is shaded gray on the x axis. Measurements of the age of the halo 
of our Galaxy yield tcai = 12.2 it 0.5 Ga, whereas measurements of the age of the disk of 
our Galaxy yield tjisk = 8-7 ± 0.4 Ga (Table 2). These age ranges are also labeled and 
shaded gray. The (Qm,f^A) = (0.3,0.7) model fits the constraints of Table 1 better than 
the other models shown. Over the past few billion years and on into the future, the rate of 
expansion of this model increases {R > 0). This acceleration means we are in a period of 
slow inflation. Other consequences of a ilA-dominated universe are discussed in (50). On 
the X axis h = 0.68 has been assumed. For other values of h, multiply the x axis ages by 
0.68/h. Redshifts are indicated on the right. 
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Fig. 2. — Age estimates of tiie universe and of the oldest objects in our Galaxy. The four 
estimates of the age of the universe from this work are indicated in Table 2. The three 
similar points near 13.4 Ga, result from h = 0.64, 0.68, 0.72 and indicate that the result 
is not strongly dependent on h when a reasonable h uncertainty of ±0.10 is used. Among 
the four, the highest value at 14.6 Ga comes from assuming h = 0.64 it 0.02. All the 
estimates in the top section of Table 2 are plotted here. As in Fig. 1, averages of the 
ages of the Galactic halo and Galactic disk are shaded gray. The absence of any single age 
estimate more than ~ 2a from the average adds plausibility to the possibly overdemocratic 
procedure of computing the variance-weighted averages. The result that to > tcai is logically 
inevitable, but the standard Einstein-deSitter model does not satisfy this requirement unless 
h < 0.55. The reference for each measurement is given under the x axis. The age of the 
sun is accurately known and is included for reference. Error bars indicate the reported la 
limits. 
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Table 1: Parameter estimates from non-CMB measurements. I refer to these as constraints. 
I use the error bars cited here as la errors in the likelihood analysis. The first four constraints 
are plotted in Fig. 3 B through E. 



Method Reference Estimate 



SNe 


{35) 


n^=o = -0.28 ±0.16 




= 0.27 ±0.14 


Cluster mass-to-light 




= 0.19 ±0.14 






Cluster abundance evolution 


en 


OA=0 _ n 1 7+0-28 
^'m — ^•J^'-O.IO 




n OO+0.25 

— U-^^_o.lo 


Double radio sources 


(8) 


qA=0 _ _n 2ti+0-70 




_ r. 1+0.50 

— ^•■'■-0.20 


Baryons 


[19) 


n^h2/3 = 0.19 ±0.12 






Hubble 


{16) 


h = 0.68 ±0.10 







Here I have combined recent independent measurements of CMB anisotropies (5), type la 
SNe {4, 5), cluster mass-to-light ratios (6), cluster abundance evolution (7), cluster 
baryonic fractions {14), deuterium-to- hydrogen ratios in quasar spectra {15), double- lobed 
radio sources {8), and the Hubble constant {16) to determine the age of the universe. The 
big picture from the analysis done here is as follows (Figs. 1 and 2): The Big Bang 
occurred at '^13.4 Ga. About 1.2 billion years (Gy) later, the halo of our Galaxy (and 
presumably the halo of other galaxies) formed. About 3.5 Gy later, the disk of our Galaxy 
(and presumably the disks of other spiral galaxies) formed. This picture agrees with what 
we know about galaxy formation. Even the recent indications of the existence of old 
galaxies at high redshift {17) fit into the time framework determined here. In this sense, 
the result is not surprising. What is new is the support given to such a young age by such 
a wide array of recent independent measurements. 



Method 



Any measurement of a function of h, fi^, and fl\ can be included in a joint likelihood 

N 



(1) 



1=1 



which I take as the product of seven of the most recent independent cosmological 

constraints (Table 1 and Fig. 3). For example, one of the Ci in Eq. 1 represents the 
constraints on h. Recent measurements can be summarized as h = 0.68 ± 0.10 {16). I 
represent these measurements in Eq. 1 by the likelihood. 



^Hubbleih) = exp 



-0.5 



h-h 
0.10 



(2) 



Another £j in Eq. 1 comes from measurements of the fraction of normal baryonic matter 
in clusters of galaxies {14) and estimates of the density of normal baryonic matter in the 
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Fig. 3. — The regions of the {Qm,^A) plane preferred by various constraints. (A) Cosmic microwave 
background, (B) SNe, (C) cluster mass-to-liglit ratios, (D) cluster abundance evolution, (E) double radio 
lobes, and (F) all combined. The power of combining CMB constraints with each of the other constraints 
(Table 1) is also shown. The elongated triangles (from upper left to lower right) in (A) are the approximate 
la, 2(7 and 3cr confidence levels of the likelihood from CMB data, £cmb {9). (A) also shows the important h 
dependence of £cmb. The contours within the dark shaded region are of h values that maximize £cmb for a 
given (fim, ^^a) pair (h = 0.70, 0.90 contours are labeled). This correlation between preferred h and preferred 
(f2m,fiA) helps £cMB(h, f2m, Oa) constrain to- In (B) through (E), thin contours enclose the la (shaded) and 
2a confidence regions from separate constraints, and thick contours indicate the la, 2a and 2>a regions of 
the combination of £cmb with these same constraints. (F) shows the region preferred by the combination of 
the separate constraints shown in (B) through (E) (thin contours) as well as the combination of (A) through 
(E) (thick contours). The best fit values are Qa = 0.65 ± 0.13 and Qm = 0.23 ± 0.08. In (A), the thin iso-to 
contours (labeled "10" through "14") indicate the age in Ga when h = 0.68 is assumed. For reference, the 
13- and 14-Ga contours are in all the panels. To give an idea of the sensitivity of the h dependence of these 
contours, the two additional dashed contours in (A) show the 13-Ga contours for h = 0.58 and h = 0.78 (the 
la limits of the principle h estimate used in this paper). In (F), it appears that the best fit has to « 14.5 Ga, 
but all constraints shown here are independent of information about h; they do not include the h dependence 

of £cMB, i^baryons Or ^Hubble (Table 1). 
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universe [Qbh^ = 0.015 ± 0.005 {15, 18)]. When combined, these measurements yield 
^m^^^^ = 0-19 i 0-12 (19), which contributes to the hkehhood through 



The ($7^, OA)-dependencies of the remaining five constraints are plotted in Fig. 3 (20). 

The 68% confidence level regions derived from CMB and SNe (Fig. 3, A and B) are nearly 
orthogonal, and the region of overlap is relatively small. Similar complementarity exists 
between the CMB and the other data sets (Figs. 3, C through E). The combination of 
them all (Fig. 3F) yields Qa = 0.65 ± 0.13 and Qm = 0.23 ± 0.08 (21). 

This complementarity is even more important (but more difficult to visualize) in 
three-dimensional parameter space: {h,Q,m,0,A). Although the CMB alone cannot tightly 
constrain any of these parameters, it does have a strong preference in the 
three-dimensional space {h,^rn,^A)- In Eq. 1, I used CcMB{h,^m-,^A), which is a 
generalization of CcMB{^m,^A) (Fig. 3A ){22). To convert the three-dimensional 
likelihood C{h, Q^, ^a) of Eq. 1 into an estimate of the age of the universe and into a 
more easily visualized two-dimensional likelihood, C{h,to), I computed the dynamic age 
corresponding to each point in the three-dimensional space (/i, fi^, ^a). For a given h and 
to, I then set >C(/t, to) equal to the maximum value of C{h, O^, Cl\) 



This has the advantage of explicitly displaying the h dependence of the to result. The 
joint likelihood C{h,to) of Eq. 4 yields an age for the universe: to = 13.4 it 1.6 Ga (Fig. 
4). This result is a billion years younger than other recent age estimates. 

What one uses for Cnubbleih) in Eq. 1 is particularly important because, in general, we 
expect the higher h values to yield younger ages. Table 2 contains results from a variety of 
h estimates, assuming various central values and various uncertainties around these 
values. The main result to = 13.4 it 1.6 Ga has used h = 0.68 =b 0.10 but does not depend 
strongly on the central value assumed for Bubble's constant (as long as this central value 
is in the most accepted range, 0.64 < h < 0.72) or on the uncertainty of h (unless this 
uncertainty is taken to be very small). Assuming an uncertainty of 0.10, age estimates 
from using h = 0.64,0.68 and 0.72 are 13.5, 13.4 and 13.3 Ga, respectively (Fig. 2). Using 
a larger uncertainty of 0.15 with the same h values does not substantially change the 
results, which are 13.4, 13.3, 13.2 Ga, respectively. For both groups, the age difference is 
only 0.2 Gy. If to oc 1/h were adhered to, this age difference would be 1.6 Gy. Outside the 
most accepted range the h dependence becomes stronger and approaches to oc 1/h (23). 

To show how each constraint contributes to the result, I convolved each constraint 
separately with Eq. 2 (Fig. 5). The result does not depend strongly on any one of the 
constraints (see "all - x" results in Table 2). For example, the age, independent of the SNe 





C{h, to) = max 



£(/t, Qrn, ^A)\t{h,nm,nA)^to 



(4) 
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data, is to{all — SNe) = 13.3l];;8 which differs neghgibly from the main result. The 
age, independent of the SNe and CMB data, is to{all - CMB - SNe) = 12.6110 Ga, 
which is somewhat lower than the main result but within the error bars. 



The Oldest Objects in Our Galaxy 

The universe cannot be younger than the oldest objects in it. Thus, estimates of the age 
of the oldest objects in our Galaxy are lower limits to the age of the universe (Table 2 and 
Fig. 2). A standard but simplified scenario for the origin of our Galaxy has a halo of 
globular clusters forming first, followed by the formation of the Galactic disk. The most 
recent measurements of the age of the oldest objects in the Galactic disk give 
tdisk = 8.7 lb 0.4 Ga (Table 2). The most recent measurements of the age of the oldest 
objects in the halo of our Galaxy give tcai = 12.2 ± 0.5 Ga (Table 2). The individual 
measurements are in good agreement with these averages. There are no large outliers. In 
contrast to the to{h, fim, Qa) estimates obtained above, all of these age estimates are 
direct in the sense that they have no dependence on a Big Bang model. 

How old was the universe when oTir Galaxy formed? If we write this as tcai + At = toi 
then what is the amount of time (At) between the formation of our Galaxy and the 
formation of the universe? If we had an estimate of At, then we would have an 
independent estimate of to to compare to to = 13.4 ±1.6 Ga, obtained above. However, we 
have very poor constraints on At. The simple but plausible estimate At 1 Gy is often 
invoked, but estimates range from ~ 0.1 to ~ 5 Gy, and may be even larger {24, 25). This 
uncertainty in At undermines the ability of estimates of the age of the oldest objects in 
our Galaxy to tell us the age of the universe. Without At, we cannot infer to from tQal- 
The best estimate of At may come from the difference between the age reported here and 
the estimate of the age of our Galaxy (Table 2) . Thus 
At = to - tGai = 13.4 - 12.2 = 1.2 ± 1.8 Gy. 

The age measurements in Table 2 also indicate that there is a 3.5-Gy period between halo 
and disk formation {toai — tdisk)- If our Milky Way is typical, then this may be true of 
other spiral galaxies. With the best fit values obtained here for the three parameters, 
(/j, n„i, Qa) = (0.72 ± 0.09, 0.23 ± 0.08, 0.65 ± 0.13), the ages tdisk and tcai can be 
converted into the redshifts at which the disk and halo formed: z^isk — l.Stli and 
ZGal = 6.0l2°3. Thus, a diskless epoch should be centered at a redshift between z^isk and 
ZGal (1-3 ~ Zdiskless ~ 6.0). We would expect fewer disks in the halolike progenitors of 
spiral galaxies in this redshift range. Studies of galaxy types in the Hubble Deep Field 
indicate that this may be the case {26). 

The requirement that the universe be older than our Galaxy, to > toah is a consistency 
test of the Big Bang model. The best fit model obtained here passes this test. There is no 
age crisis. This is true even if the high values of h {^-^ 0.80) arc correct. Only at h ~ 0.85 
is to ~ tcai- This consistency provides further support for the Big Bang model, which the 
standard model ($7^ = 1, Oa = 0) is unable to match unless h < 0.55. 
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Fig. 4. — This plot shows the region of the h — to plane preferred by the combination 
of all seven constraints. The result, to = 13.4 it 1.6 Ga, is the main result of this 
paper. The thick contours around the best fit (indicated by a star) are at likelihood levels 
defined by £/£max = 0.607, and 0.135, which approximate 68% and 95% confidence levels, 
respectively. These contours can be projected onto the to axis to yield the age result. This 
age result is robust to variations in the Hubble constraint as indicated in Table 2. The areas 
marked "Excluded" (here and in Fig. 5) result from the range of parameters considered: 
0.1 < fim < 1-0, < Qa < 0.9 with + < 1- Thus, the upper (high to) boundary is 
defined by (Qm,fiA) = (0.1,0.9), and the lower boundary is the standard Einstein-deSitter 
model defined by (Qm,f^A) = (1)0). Both of these boundary models arc plotted in Fig. 1. 
The estimates from Table 2 of the age of our Galactic halo (tcai) and the age of the Milky 
Way (tdisk) are shaded grey. The universe is about 1 billion years older than our Galactic 
halo. The combined constraints also yield a best fit value of the Hubble constant which can 
be read off of the x axis (h = 0.73 it 0.09, a slightly higher and tighter estimate than the 
input h = 0.68 ± 0.10). 
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Fig. 5. — The purpose of this figure is to show how Fig. 4 is built up from the seven 
independent constraints used in the analysis. All six panels are analogous to Fig. 4 but 
contain only the Hubble constraint [h = 0.68 ± 0.10, (Eq. 2)] convolved with a single 
constraint: (A) cosmic microwave background, (B) SNe, (C) cluster mass-to-light ratios, 
(D) cluster abundance evolution, (E) double radio lobes, and (F) baryons (Table 1). The 
relative position of the best fit (indicated by a star) and the 13.4-Ga line indicates how each 
constraint contributes to the result. 
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Comparison with Previous Work 

The goal of this paper is to determine the absolute age of the universe toih, ^a)- 
Knowledge of h alone cannot be used to determine to with much accuracy. For example, 
the estimate h = 0.68 ± 0.10 corresponds to 8 < < 22 Ga (Fig. 4). Similary, knowledge 
of (Q^m^A) yields Hoto{^m,^A) not to {Ho is the usual Hubble constant). When one 
inserts a preferred value of h into a Hoto result, one is not taking into consideration the 
correlations between preferred h values and preferred {^lm,^A) values that are inherent, 
for example, in CcMB{h, O^, and Cbaryonsih, ri^)- The preferred values of h in these 
likelihoods depend on and Q,a- Perlmutter et al.{4) used SNe measurements to 
constrain {Clrn,^A) and obtained values for Hotg- To obtain to, they did the analysis with 
h set equal to the value preferred by their SNe data, h = 0.63. Their result is 
to = 14.5 lb 1.0(0.63//i) Ga. When a flat universe is assumed, they obtain 
^flat ^ i4.9+i-4(o.63//i) Ga. Riess et al.{5) found h = 0.65 ± 0.02 from their SNe data. 
Marginalizing over this Hubble value and over JIa and Q.ra, they report to = 14.2 ±1.7 Ga. 
When a flat universe is assumed, their results yield t^'"* = 15.2 it 1.7 Ga. The Perlmutter 
et al.{4) and Riess (5) results are in good agreement. When I assume h = 0.64 it 0.02, I 
get to = 14:.6~^\'i Ga. This result is plotted in Fig. 2 to illustrate the important influence 
on the result of using a small h uncertainty. Efstathiou et aL(12), on the basis of a 
combination of CMB and Perlmutter et al.(4) SNe data, have estimated 
to = 14.6(/i/0.65)~-'^ Ga. I used h = 0.65 ± 0.0 with this data combination to get 
to = 14.5^ j^lo Ga. However, when I used h = 0.65 ± 0.10, the result is 0.7 Gy lower 
{to = 13.8lf j Ga). To obtain the main result, I used uncertainties large enough to reflect 
our knowledge of h on the basis of many sources. The use of a larger h uncertainty 
contributes to the substantially younger ages found here {23). 

A potential problem with the SNe ages is the high region, {0,^, ^a) ^ (0-8, 1.5), which 
dominates the SNe fit. This region is strongly disfavored by the six other constraints 
considered here (see Fig. 3). These high {Q,jn,^A) values allow lower ages than the i^'"* 
SNe results because the slope of the iso-to contours (Fig. 3B) is larger than the slope of 
the SNe contours. The t^'"* results are not as subject to this problem and are the results 
most analogous to the result reported here, despite the fact that the SNe i^'"* results are 
less consistent with the result reported here. There are several independent cosmological 
measurements which have not been included in this analysis cither because a consensus 
has not yet been reached [gravitational lensing limits {27-30)] or because the analysis of 
the measurements has not been done in a way that is sufficiently free of conditioning on 
certain parameters [local velocity field limits {31)]. Doubts about some of the observations 
used here are discussed in {32). There has been speculation recently that the evidence for 
r^A is really evidence for some form of stranger dark energy that we have incorrectly been 
interpreting as Several workers have tested this idea. The evidence so far indicates 
that the cosmological constant interpretation fits the data as well as or better than an 
explanation based on more mysterious dark energy {4, 33, 34). 
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Table 2: Age estimates of our Galaxy and universe (36). "Technique" refers to the method used to make 
the age estimate. OC, open clusters; WD, white dwarfs; LF, luminosity function; GC, globular clusters; M/L, 
mass-to-light ratio; and cl evol, cluster abundance evolution. The averages are inverse variance-weighted 
averages of the individual measurements. The sun is not included in the disk average. "Isotopes" refers 
to the use of relative isotopic abundances of long-lived species as indicated by absorption lines in spectra 
of old disk stars. The "stellar ages" technique uses main sequence fitting and the new Hipparcos subdwarf 
calibration. "All" means that all six constraints in Table 1 and the CMB constraints were used in Eq. 1. 
"AU-x" means that all seven constraints except constraint x were used in Eq. 1. Figures. 3 and 5 and the all 
- x results indicate a high level of agreement between constraints and the lack of dependence on any single 
constraint. Thus, there is a broad consistency between the ages preferred by the CMB and the six other 
independent constraints. Figure 2 presents all of the disk and halo age estimates. 



Technique 


Reference 


h Assumptions 


Aee (Ga) 

' o v. / 


Object 


Isotopes 




None 


4.53 ± 0.04 


Sun 


Stpllar affps 


(38) 


None 


8.0 ± 0.5 


Disk OC 


WD LF 


(39) 


None 


8.0 ± 1.5 


Disk WD 




(in) 


None 


9.0 ± 1 


Disk OC 


WD LF 


(9^) 


None 


Q 7+0.9 
^■'-0.8 


Disk WD 


Stpllar affps 


(Al) 


None 


12 O+o S 


Disk OC 






None 


8.7 ± 0.4 


X A- 1 ( TWfPl 




(J,9) 




11.5 ± 1.3 


Halo GC 


^kdicii agC9 


(J1) 


None 


11 8+^-^ 


Halo GC 

1 ICI vJ Vw. 


^LCiicii agC9 


(J J) 


None 


12 ± 1 


Halo GC 


Stellar ages 


(45) 


None 


12 ±1 


Halo GC 


Stellar ages 


m 


None 


12.5 ± 1.5 


Halo GC 


Isotopes 


m 


None 


13.0 ±5 


Halo stars 


Stellar ages 


m 


None 


13.5 ± 2 


Halo GC 


Stellar ages 




None 


14.0li:^ 


Halo GC 






None 


12.2 ±0.5 


tcai (avg) 


SNe 


U) 


0.63 ± 0.0 


14.5 ± 1.0 


Universe 


SNe (flat) 


U) 


0.63 ±0.0 




Universe 


SNe 


(5) 


0.65 ± 0.02 


14.2 ± 1.7 


Universe 


SNe (flat) 


(5) 


0.65 ± 0.02 


15.2 ± 1.7* 


Universe 


All 


This work 


0.60 ±0.10 


15.51-;: 


Universe 


All 


This work 


0.64 ± 0.10 


13.512-2* 


Universe 


All 


This work 


0.68 ±0.10 


13.4l}:«* 


Universe 


All 


This work 


0.72 ±0.10 


13.3l}f,* 


Universe 


All 


This work 


0.76 ± 0.10 


12.31- 


Universe 


All 


This work 


0.80 ±0.10 


11.91}:^ 


Universe 


All 


This work 


0.64 ± 0.02 


14.6+}-?* 


Universe 


All-CMB 


This work 


0.68 ± 0.10 


14.01^:" 


Universe 


AII-SNe 


This work 


0.68 ±0.10 




Universe 


All-M/L 


This work 


0.68 ±0.10 


13.3+}? 


Universe 


All-cl evol 


This work 


0.68 ± 0.10 


13.3l}I 


Universe 


All-radio 


This work 


0.68 ±0.10 


13.31}:! 


Universe 


All — baryons 


This work 


0.68 ±0.10 


i3.4i?:« 


Universe 


All-Hubble 


This work 


None 


< 14.2 


Universe 


AII-CMB-SNe 


This work 


0.68 ± 0.10 


12.61^0 


Universe 



* Also plotted in Fig. 2. 
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